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Abstract. - 

We study the effect of magnetic scattering on transport in a system with strong structural 
disorder, using exact finite size calculation of the low frequency optical conductivity. At weak 
electron-spin coupling, spin disorder leads to a decrease in resistivity, by weakening the quantum 
interference precursors to Anderson localisation. However, at strong electron-spin coupling, the 
'double exchange' limit, magnetic scattering increases the effective disorder, sharply increasing 
the resistivity. We illustrate the several unusual transport regimes in this strong disorder 
problem, identify a re-entrant insulator-metal-insulator transition, and map out the phase 
diagram at a generic electron density. 



Introduction. - The physics of transport and localisation in the presence of structural 
disorder has been extensively studied [1]. In three dimension (3d) increasing disorder leads 
to a monotonic increase in the fraction of localised states, and the resistivity in the extended 
regime, leading finally to the Anderson metal- insulator transition (MIT). 

The interplay of scattering from 'paramagnetic' moments with scattering from structural 
disorder leads to several new transport regimes, whose character is poorly understood. The 
presence of weak magnetic scattering actually weakens the localising effect of disorder, as 
discovered by Lee [2] and by Hikami et al. [3], while strong magnetic coupling, the 'double 
exchange' limit, enhances the localising effect of structural disorder [4] . There is no understand- 
ing of how these two endpoints are connected. This experimentally relevant "middle" is wide, 
unexplored, and beyond the reach of standard Boltzmann transport theory [5] . In this paper we 
present essentially exact results on the conductivity and MIT considering the combined effect 
of structural disorder and magnetic scattering, and illustrate the novel transport regimes in 
the problem. 

An understanding of transport properties of disordered magnetic systems is of particular 
relevance now because of intense experimental activity in diluted magnetic semiconductors 
[6, 7, 8] (DMS), amorphous magnetic semiconductors [9, 10], e.g GdSi, and the manganites 
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[11]. Some of these systems, notably those where the magnetic coupling arises from Hunds 
rule, as in d electron systems, require an understanding beyond the 'weak magnetic scattering' 
studies in the localisation literature. A quick look at the resistivity in the paramagnetic phase 
in these materials reveal that they are all 'poor metals'. The resistivity at 300K in the DMS, 
Ga 1 _ x Mn a; As, at x <~ 0.08, is ss 4 — 6 mOcm [8]. For a-GdSi, in its 'metallic regime' this is 
<~ 3 intern [9]. For the manganite Lai-^Sr^MnOs (LSMO), at large doping, x > 0.3, where 
electron-phonon coupling effects are expected to be weak, the room temperature resistivity is 
> 5 mOcm [12]. To put these numbers in perspective, let us use the Mott 'minimum metallic 
conductivity' [13] as reference. The Mott 'minimum' conductivity is ~ 0.03(e 2 /^ a o) and if we 
use a lattice constant do = 3A then UMott ~ 2.5 x 10 4 (Qm) _1 , i.e, the 'Mott resistivity' pMott 
is approximately 4 mficm. Comparing with the experimental data quoted above, for all the 
systems concerned, p/pMott ~ 0(1) in the paramagnetic phase. There is no standard theory 
for analysing the scattering from structural and magnetic disorder in this regime. 

A complete theory of transport, for example the temperature dependence, in any of these 
systems would of course have to deal with annealed spin disorder. We will discuss these effects 
in the future, but focus here on the simpler case of scattering from a combination of quenched 
structural and magnetic disorder. Our principal results are contained in Fig.l, showing the 
'global' behaviour of the resistivity, and the phase diagram in Fig. 3. 

Model. - We study transport in the following model: 

H = -t^2 c\ a c jt7 + ^2(ci - /j,)n i(7 - J'^ai.Si (1) 

{ij),a ia i 

This is the 'Anderson disorder' problem, with the electrons being coupled additionally to 
spins, Sj, through J'. We use J' > 0. The hopping is only between nearest neighbour sites in a 
simple cubic lattice. The random on site potential, ej, is uniformly distributed between ±A/2. 
The dimcnsionless parameters in the problem are disorder A/t, magnetic coupling J'S/t and 
the electron density n, controlled by p. We absorb S in our magnetic coupling J', assuming 
|Sj| = 1. W — I2t is the bandwidth. We assume that the spin distribution is uncorrelated 
between sites, and isotropic on site, i.e each spin can point anywhere on the surface of a unit 
sphere. The energies are measured in units of t, set to 1. 

Conductivity calculation. - We estimate the d.c conductivity, <Jdc, by using the Kubo- 
Greenwood expression [14] for the optical conductivity. In a disordered non interacting system: 

ff H = 4 £("" - np)^Ls(w (ep - e«)) (2) 

The constant A — (ire 2 )/hao. The matrix element f a p — (4>a\jx\ipp) and the current operator 
in the tight binding model is j x = itaoej^i <j( c l+xa <r c i,<? ~ h.c). The tp a etc arc single particle 
eigenstates, for a given realisation of disorder, and e a , ep etc are the corresponding eigenvalues. 
The n a = 9(p — e a ) etc are occupation factors. 

The conductivity above is prior to disorder averaging. We work in a Lt x Lt x L geometry 
and, given the finite size, the 5 function constraint in a(uj) cannot be satisfied for arbitrary oj. 
However, we can still calculate the average 'low frequency' conductivity, defined below, and 
extract the d.c. conductivity by studying a sequence of system sizes. 

We set Lt = 6, and sum over the 5 functions to compute the low frequency average 
cr av (p, Alu, L) = (Alu)^ 1 J o a " a(p, uj, L)duj, for {L : 24, 28, 32, 36, 40}. using periodic boundary 
condition in all directions. The averaging interval is reduced with increasing L, with Auj ~ 
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Fig. 1. - Variation of resistivity, p, with J' and A. Density n = 0.5. Data obtained by extrapolation 
on a sequence of sizes 6 x 6 x L, with L varying from 24 to 40. Panel (a) shows p in the intermediate 
J' regime, while panel (b) shows a 'global' view, with a logarithmic J' scale. Panel (c) shows how the 
'd.c conductivity' is obtained by extrapolation, illustrating the behaviour for various J' at A = 16, 
i.e, near the Anderson MIT. The conductivity in panel (c) is normalised to a(L = 24). 

B/L. The constant B is fixed by setting Alu = 0.06 at L = 40. a av is disorder averaged 
over N r realisations, which involve 'cross averaging' over ej and spin configurations, i.e both 
the ej and Sj are chosen freshly for each realisation. We use N r ~ 100 for the largest size 
(L = 40), increased to N r — 400 for L = 24. Denoting the disorder averaged low frequency 
conductivity for a size L and Auj — B j ' L, as o~ av (/j,, B/L, L), the d.c conductivity is obtained as 
Vdc{v) = Hm.L^ooO'avi^, B/L, L). The chemical potential is set to target the required electron 
density n. Our transport calculation method and some benchmarks will be discussed in detail 
elsewhere [15]. To convert to 'real' units, note that our conductivity results are in units of 
(ne 2 )/hao, i.e a = 1 on our scale corresponds to w 10 2 <JMott- All our results in this paper are 
for n = 0.5, which is 'quarter filling' in the weak J' limit, and half-filling (of the lower band) 
at large J', n = 0.5 is a generic 'high density' case and exhibits all the interesting regimes in 
transport. The density dependence of transport properties will be discussed separately [15]. 

Transport Regimes and Insulator- Metal Transitions. - The qualitative features in transport 
are immediately visible in the resistivity, Fig.l, panels (a) — (6), the most noteworthy being 
the strongly non monotonic behaviour of p(J') at large A. To set the stage for discussing the 
interplay of structural disorder and magnetic scattering, let us quickly review the individual 
effects arising from each of them. 

Analytic approaches to the structural disorder problem, beyond Born scattering, consider 
the weak localisation (Cooperon) contributions to the conductivity in a (fc^/) -1 expansion, 
with Uf and I being the Fermi wavevector and mean free path respectively. Beyond this 
one can use the approximate self-consistent theory (SCT) of localisation [16], or depend on 
numerical results. With growing disorder the resistivity, at fixed electron density, increases 
monotonically, finally diverging at the Anderson transition. 

For magnetic disorder, the lowest order effect arises from Born scattering, with p(J') cx J' 2 . 
Beyond this regime the resistivity shows a dependence [15] of the form p(n, J') ~ bi(n)J' 2 + 
62 (n) J' 4 , upto moderate coupling, J'/t ~ 3, with 61,62 being electron density dependent 
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Fig. 2. - Left Fig: Resistivity at small J': (a) Quadratic fit to the extrapolated resistivity: p(J') ~ 
pi (A) + ci(A)| J'| + di(A)J' 2 , at weak coupling, and moderate A. (b) pi (A), with a polynomial fit 
(see text), (c) The coefficient — ci, with a power law fit (see text), and (d) The coefficient d\. Right 
Fig: Resistivity at large J': (e) Fit to a(J' , A) of the form a(oo, A) + fi/J'- (/) Resistivity in the 
double exchange limit p(oo, A) = l/cr(oo, A), (g) The coefficient f\. 



coefficients, At strong coupling, and the band center, p(J') saturates with growing J'. 

Our interest is in the combined effect of these two scattering mechanisms. There are 
tentatively five different transport regimes in the problem, and, except for (i) below, all of 
them are beyond the reach of standard transport theory. These are: (i) the weak scattering 
regime, where Mathiessens rule holds, (ii) spin flip correction to weak localisation (WL), with 
the A dependence showing WL corrections and spin flip increasing the conductivity, (Hi) spin 
dephasing driven insulator-metal transition (IMT), occuring over a window in A, (iv) the 
disordered double exchange (DE) limit, and (v) the intermediate coupling 'metal'. We analyse 
our data in the spirit of this classification. 

(i) Weak scattering, obeying Mathiessens rule: The regime of weak structural disorder and 
weak electron-spin coupling can be understood in terms of additive Born scattering, with the 
net scattering rate, T(A, J') w aiA 2 + 6i J' 2 . The resistivity is additive and both contributions 
are described by lowest order perturbation theory This corresponds to the bottom left hand 
corner in Fig.l.(a). The 'parallel' curves in Fig. 2 (a) for A upto 4, and the regime J' < 3 
broadly identify the domain of 'Mathiessens rule'. The net resistivity in this regime, Fig. 2. (a), 
is < O.lpMott, i-e, a few hundred p;£!cm. 

(ii) Spin flip correction to weak localisation: In the structural disorder problem the quantum 
corrections to the conductivity become important with growing A, and show up via the WL 
effect. This effect is already beyond standard transport theory. The leading corrections to 
cr(A) in three dimension, beyond Boltzmann transport, have been worked out [17]. Exact 
numerical calculations [15, 18] suggest that the Boltzmann result and the low order quantum 
corrections describe the resistivity upto A/W w 1. In this regime, with A/W ~ 1, the effect 
of magnetic scattering, with J' /W <C 1 is non trivial. Spin flip scattering of the electrons 
[2, 3], by the random magnetic moments serves to suppress the localising effect of structural 
disorder, Fig.l.(a) and Fig. 2. (a), as we discuss below. 

Just as inelastic scattering leads to decoherence, and a suppression of quantum interference, 
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spin flip scattering of the electrons, arising from processes of the form (S xi + iS^c^Cjf etc, 
leads to spin decoherence and a corresponding cutoff to WL corrections. If the dephasing time 
is T a , then the spin diffusion length is l s = y/Dr s , with the diffusion constant D oc A~ 2 , and 
t s -1 oc J' 2 . The WL correction gets corrected [1] by da oc I" 1 oc A|J|. If the WL correction to 
the conductivity scales as <~ (/cfO 1 ~ A 2 , then the corresponding spin flip correction to the 
resistivity should behave as Sp ~ — A 5 | J'\. 

Fig. 2. (a) shows p(J',A) in the moderate J' region, with A increasing beyond the Born 
scattering regime. At weak disorder the resistivity has the expected form p(A,J') oc aiA 2 + 
&i J' 2 , while at larger disorder we expect the A dependence to be stronger, and simultaneously 
a positive contribution to the conductivity 8a oc \J'\ to show up. The coefficient of the J' 2 
term could also be renormalised at large A. Taking these possibilities into account we fit the 
low J' resistivity to the form: p(A, J')| j'-o ~ Pi (A) + h J' 2 + Ci(A)| J'\ + /(A) J' 2 , where the 
coefficients would depend on electron density. p\ (A) tracks the localising effect of structural 
disorder, ci(A) measures the 'antilocalising' effect of spin flip scattering, and d\ = b\ + /(A) 
monitors the renormalised 'Born scattering' from spins. 

For n = 0.5, the resistivity from structural disorder, pi(A) w aiA 2 + a 2 A 4 + a 3 A 6 , as we 
fit in Fig.2.(b), with a x = 0.408, a 2 = 3 x 10~ 4 , and a 3 = 4 x 10~ 5 . The coefficient of the 
'anti- localisation' term, Fig.2.(c) is ci(A) w coA", with co = —3.4 x 10~ 4 , and a <~ 5. The 
(renormalised) coefficient of the J' 2 term, di = b\ +/(A), has the form indicated in Fig.2.(d). 
While Fig.2.(b) highlights the rapid growth in resistivity due to structural disorder, Fig.2.(c) 
indicates how this growth is cutoff via the spin flip scattering effected by J'. This affects the 
metal- insulator phase boundary as we discuss next. 

(Hi) Spin dephasing driven IMT: There is no perturbative scheme as the disorder increases 
beyond the weak localisation regime. The SCT provides a guide and of course there is extensive 
numerical work on the Anderson transition. Unfortunately, the SCT has no equivalent that 
includes magnetic scattering as well (although there were early attempts in two dimension [19]), 
and we do not know of even numerical work probing the MIT including magnetic scattering. 

As A drives the system towards localisation, the effect of weak J' is dramatic. The small 
increase in conductivity, correcting the weak localisation contribution, now develops into a full 
fledged anti- localisation effect, enlarging the domain of the metallic phase, see phase diagram in 
Fig. 3. (a), and pushing up the critical disorder needed for localisation. The critical disorder now 
depends on J' and increases from A w 16.5, at J' = 0, to A <~ 27 at J' ~ 8, before dropping 
to A ~ 12 for J' — > oo (see Fig. 3. (a)). The optical conductivity data, Fig.3.(b) illustrates 
how the low frequency conductivity increases as J' grows from 0.2 to 2.0, at A = 16, and 
drastically reduces again for J' = 20. There is an insulator-metal-insulator transition with 
increasing J' for 26 > A > 16.5. 

While the IMT driven by 'small J" may be motivated as the 'extrapolation' of the antilocal- 
ising effect of spin flip scattering, new physical effects emerge rapidly as J' is increased, staying 
in the large A regime. The reducing resistivity, Fig.l.(a), and Fig. 2. (a), exhibits a minima, 
and then rises again with increasing J'. Such behaviour can be viewed as an extension of the 
J' 2 term seen at weak disorder, but it is probably more fruitful to approach the problem from 
the strong coupling, double-exchange, end. 

(iv) Double exchange with disorder: This regime corresponds to J'/W — > oo, in the presence 
of arbitrary structural disorder, i.e, the right edge of panel l.(b). The form of the resistivity 
p(J') arising from 'magnetic disorder' at large J' is very different from what one observes in 
p(A) at large A. This is because J' contributes to both 'band splitting' and effective disorder, 
and the effective disorder saturates as J'/W oo with J' controlling only the band splitting. 

A standard transformation [20] and projection maps on the t — J' problem to a spinless 
fermion model with hopping amplitudes Uj dependent on nearest neighbour spin orientation: 
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Fig. 3. - (a). The insulator-metal phase diagram in A — J', (6) the optical conductivity illustrating 
the 'I-M-F feature, (c) Combining the weak coupling and strong coupling expansions, Fig.2.(a) and 
Fig.2.(e), for a 'global' description of the resistivity. 



H = Eij Uj (9, <t>)ihj + J2 t tiihi = Eij toTiij + Ey tUjihj + E 4 tula with the hopping 

amplitude specified by tij = —t(cos^cos^ + sin^sin^-e 1 ^* - ^^). We have split this into the 
'mean hopping', to — (l^jl)i an d the fluctuation St. 

The localisation properties of this model have been studied by Li et al. [4], although they 
did not calculate the resistivity. The 'hopping disorder' by itself localises less than 0.5% of the 
states in the band, and, as we observe, the resistivity at band center remains finite, ~ 0.2pMott- 
On adding structural disorder the mobility edge moves inward, with localisation of the full band 
occuring at A/t <~ 11.5. The critical disorder, visible approximately at large J' in Fig. 3. (a), 
can be motivated by the band narrowing in to compared to t. In the fully spin disordered 
phase to/t <~ 2/3, and a crude scaling of A C (J' = 0) would suggest, A c (oo)/A c (0) <~ 2/3, so 
that the double exchange model in the paramagnetic phase would localise at A ~ 11. The 
small deviation from this value arises from the presence of 5Uj and the Berry phase involved 
in it. Apart from the mobility edge calculation [4] we do not know of any exact results on 
transport in the strongly disordered double exchange model. 

The growth in resistivity with A, remaining at large J', is visible in the right edge of 
Fig.l.(b). Within our calculation the MIT at large J' occurs slightly above A = 12. 

(v) Strongly disordered, intermediate coupling metal: This is the most complicated regime 
in the problem and also the most relevant for real disordered magnetic systems. The regime 
corresponds to A/W ~ 1 and J' /W <~ 1, and, as obvious from Fig.l, p/pMott > 1- This 
is a nominally 'metallic' but diffusive, highly resistive regime. There are no analytical tools 
for directly estimating transport properties in this regime, but much of the physics can be 
motivated by an expansion about the double exchange limit. 

When J' is finite, both the electron spin states, polarised parallel and anti-parallel to the 
core spin, need to be retained. In terms of these states the Hamiltonian assumes the form: 
H = E»j ^llaliP + E» £ i n » ~ 4" Hi( n ia — n if))- At large J' the chemical potential will be in 
the lower band, if we want n < 1. We can split the tf? , as we did for the single band model, 
into mean amplitudes and fluctuations. The major source of disorder is still ej, with additional 
contribution from the 5t"? . The orbital mixing effect of 'off diagonal' couplings, either in terms 
of mean amplitude or fluctuations, is regulated by the large energy denominator J'. Although 
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the 'reference' problem, J' — > oo, is not analytically tractable in the presence of structural 
disorder, we expect orbital mixing to generate corrections to conductivity ~ 0(1/ J'). Fig.2.(e) 
fits the data, for A upto 8 and J' down to 6, to the form a(J', A) | j'^oo ~ er(oo, A) + d\/J'. 
Fig.2.(f) show p(oo, A) = l/er(oo, A) and Fig.2.(g) is the coefficient /i. 

Fig. 3. (a) puts together the phase diagram that emerges from tracking the conductivity 
across the full A — J' range. While the A dependence seems intuitive, in that increasing 
structural disorder always leads to a MIT, the effect of J', at fixed large A, is nontrivial. We 
can identify a clear FM-I transition which, we believe, is a new result. As discussed before, 
the <j(u>) in panel 3(b) confirms the I-M-I transition seen in the resistivity. 

Fig.3.(c) puts together the weak coupling and strong coupling fits to p(J',A), for A < 8, 
comparing with the full numerical data. The quality of the fit indicates that the dominant 
transport regimes in this problem can be classified as we had suggested earlier and some of 
the physics can be motivated by extending concepts that are already available in the transport 
literature. The regions of MIT cannot be described by these expansions and should probably 
be accessed via a complementary scaling approach. 

Conclusions. - We have presented exact results on electron transport in the background 
of arbitrary structural and spin disorder and provided a framework within which the data 
can be analysed. We have identified the distinct transport regimes in the model, explored the 
Anderson transition in the presence of weak spin disorder, and located a novel insulator-metal- 
insulator transition driven by increasing magnetic coupling. In addition to the intrinsic interest 
as a non trivial disorder problem, our results should be useful in understanding transport in 
the paramagnetic phase in several currently interesting magnetic materials. 
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